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We develop a new approach to characterizing the morphology of rough surfaces based on the 
analysis of the scaling properties of contour loops, i.e. loops of constant height. Given a height 
profile of the surface we perform independent measurements of the fractal dimension of contour 
loops, and the exponent that characterizes their size distribution. Scaling formulas are derived and 
used to relate these two geometrical exponents to the roughness exponent of a self-affine surface, 
thus providing independent measurements of this important quantity. Furthermore, we define the 
scale dependent curvature and demonstrate that by measuring its third moment departures of the 
height fiuctuations from Gaussian behavior can be ascertained. These nonlinear measures are used 
C/5 ' to characterize the morphology of computer generated Gaussian rough surfaces, surfaces obtained 

in numerical simulations of a simple growth model, and surfaces observed by scanning-tunneling- 
microscopes. For experimentally realized surfaces the self-affine scaling is cut off by a correlation 
length, and we generalize our theory of contour loops to take this into account. 
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I. INTRODUCTION 



o. 

O , Random surfaces are widely used in the physical sciences to model phenomena ranging from the extremely small 
(quantum gravity) to the very large (Earth's relief). They describe crack fronts in materials science ripple- wave 
turbulence passive tracers in two-dimensional fluid flows (|,^, cloud perimeters shapes of stromatolites 

^ (conjectured fossil accretions of ancient bacteria) to mention but a few recent examples. Our focus in this paper 
0^ is on the morphology of deposited metal films, which develop random self-affine surfaces under several quite different 
non-equilibrium growth conditions, as indicated by theoretical, numerical, and experimental results over the past 
decade 0Jll[. 

Surface configurations are parametrized by a two-dimensional field ft.(x) which represents the height of the surface 
above a reference plane {x}. Theoretically the dynamics of a growing surface are described by a continuum (Langcvin) 
" equation giving dh{:x.)/dt as a sum of a Gaussian white noise term, to mimic the random deposition of atoms, and a 
"Y^i polynomial of various gradients of /i(x), to model relaxation processes on a coarse grained scale. The nonequilibrium 
' growth behavior is due to the interplay of the deposition and relaxation terms. 

, Different relaxation terms are appropriate for different growth conditions, and the notion of universality has been 
I ' taken over from the study of critical phenomena. Namely, it is believed that there are only a few distinct universality 
' O [ classes of growth characterized by exponents which describe the temporal and spatial scaling of the growing interface. 

One would guess that even a snapshot of the morphology should carry the evidence of the non-equilibrium, nonlinear 
growth process which produced it, and should differ measurably from surfaces produced in equilibrium or by a linear 
process, even if they share the same scaling exponents. 
' This motivates a search for roughness measures independent of the quadratic ones (e.g. the rms height), which 
, might identify important distinctions between different surface models that have similar spatial power spectra. Such 
^ measures, although motivated in the context of self-affine or multi-affine surfaces, should be handy even for surfaces 
5^ , showing no self-affine regime. They can quantify features of morphology which are presently characterized by eye, 
which should permit a more systematic comparison between observations and models than at present. One can 
imagine that, armed with two or three kinds of roughness measures tuned to different qualitative aspects of the 
surfaces morphology, one could construct empirical "phase diagrams" in this two or three-dimensional parameter 
space, e.g. mapping out domains in the parameter space that correspond to various growth conditions. 

Given the surface as parametrized by an array of heights - obtained, e.g., from a simulation or a scanning-tunneling- 
microscope (STM) experiment - we ask, in what different ways can the surface morphology be characterized? In 
general, one requires more than one characterization to confirm a match between experimental and simulation data, 
or to convincingly verify self-affineness. For the applied problems of growing fiat surfaces (e.g. for semiconductor 
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devices) or regularly modulated ones (e.g. to nanofabricate arrays of quantum dots), it is also desirable to develop 
independent measures that quantify different aspects of a rough surface's geometry. In this paper, we propose two 
categories of novel measure for characterizing spatial correlations of rough surfaces. These measures are usable on any 
kind of rough model, and require no dynamical information, so they should be useful in analyzing not only experiments 
on solid films, but all the diverse phenomena mentioned above. 



A. Outline of the paper 



We start with a short review (Section ||) of sclf-afRne geometry in terms of its real-space, Fourier space and fractal 
properties. In section [II we introduce the non-linear measures, the scale dependent curvature and the loop measures. 



for which various scaling relations are derived in section IV, These scaling relations are modified in the presence of a 
(time dependent) cutoff length-scale above which the height fluctuations are no longer self-affine, and this situation 
is described towards the end of section IV. This concludes the first half of the paper which deals with the theory of 
non-linear measures. 

The second half of the paper is devoted to analyzing data obtained in numerical simulations and experiments using 
the measures introduced in the first part. It starts off with section ^ where we present the results of our simulations of 
random Gaussian surfaces for various values of the roughness exponent < a < 1. These simulations serve to confirm 
the various scaling relations derived earlier. In section ^ the non-linear measures are applied to a non-equilibrium 
growth model, the so-called single step model, which is known to belong to the KPZ universality class. Finally, in 
section VII we demonstrate the usefulness of our measures for analyzing experimental data on the example provided 
by an STM image of a growth roughened metal film. The discussion section (section VII] ) summarizes our main 
results, gives a critical comparison between the newly introduces measures and those used previously, and points out 
some interesting new directions in which progress can be made. The three appendices are reserved for details of the 
calculation of the loop correlation exponent in the case of equilibrium rough surfaces (appendix details of the 
derivation of percolation exponents for contours of uncorrelated heights (appendix B|), as well as a full description 
of the loop finding algorithm which is at the heart of the numerical simulations and the loop analysis of STM data 
(appendix ^). 



II. SELF-AFFINE GEOMETRY 



Here we review the scaling properties of self-affine interfaces in real-space and in Fourier-space, as well as the fractal 
geometry of their level sets. The surface is fully specified by the height field /i(x), which may be the microscopic 
heights of individual surface atoms above the substrate, as measured by an STM, or it may be a coarse grained 
quantity representing the average of individual atomic heights over a region fl^ . 

The defining property of self-affine surfaces is their invariance under rescaling. Namely, the probability distribution 
function (PDF) for h{x.) is such that 

/i(x) = b-°'h{bx) (2.1) 

for any 6 > 1, where a is the roughness exponent; here the symbol = means "statistically equivalent with respect to 
the PDF". In other words, if we stretch the surface by a rescale factor b in the horizontal direction (= parallel to 
the reference plane x), then to obtain a statistically equivalent surface, we must stretch by a factor 6" in the vertical 
direction (= the perpendicular direction of the heights h{x.)). A central theme of this paper is the different ways of 
determining a, given a height profile /i(x). 

A self-affine surface is rough if a > 0. Furthermore, for the surface to exhibit a two-dimensional character (at 
distances much larger than the surface width) the rescale factor in the vertical direction {b") can not exceed the one 
in the horizontal direction (6), i.e., we require a < 1. 



A. Real-space properties 



The self-affine scaling of the height is typically measured by the height-correlation function 

i?2(r) = ([Mx + r)-Mx)]2)~|rp" , 



(2.2) 



where the scaling with separation |r| is a direct consequence of the self-affine property, Eq. (2.1), which states that h 
has a scaling dimension a. 
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In experiments the correlations that lead to self-afRne scaling of the surface develop over time, which we take to be 
measured from the start of the deposition process. Namely, after time t self-afRne scaling will be observed only up-to 
length scales smaller then the correlation length ^{t). Physically the height correlations develop due to the various 
surface relaxation processes that are present under the given growth conditions. 

Numerical simulations of various surface growth models, as well as experiments under different conditions, have 
shown that ^ grows with the duration of the deposition process, t, according to the dynamical scaling relation |p^ : 

m - t'^' ■ (2.3) 

It is believed that there are only a few different universality classes of growth each characterized by the exponents a 
and z ^1 14 1 . Experimental efforts have been focused on extracting these exponents from data obtained using various 



surface-sensitive methods: X-rays or helium diffraction, STM scans, etc. |10| In this paper, we will be almost entirely 
concerned with the spatial (equal-time) correlations. 

B. Fourier-space properties 

The power spectrum of a self-afRne surface 

^(q)-(IMq)l') (2.4) 
is defined in terms of the Fourier transformed height 

h{^) ^ J H^K'"^"" ■ (2.5) 

The height correlation function is linearly related to the height power spectrum, 

D2{r)^ fd\S{q){e'^''-l), (2.6) 



as is any other translation-invariant expectation quadratic in heights, such as the net variance. Eq. (2^) and Eq. (2.2) 
imply the scaling 

5(q) ~ |q|-2(i+") , (2.7) 
for small values of |q|. In the case of a surface with a finite corr elatio n length, S'(q) crosses over to a constant value 



for |q| < 1/C(i); this situation is discussed in more detail in Sec. IV E , 

Clearly, S{q} does not uniquely characterize a self-afRne ensemble of surfaces. For example, it is invariant under 
/i(x) — !■ — /i(x), yet surfaces produced in non-equilibrium growth typically break the up/down symmetry. Furthermore, 
given any S'(q) one can always construct a Gaussian ensemble by linear addition of Fourier components - we do this 
in Sec. ^ - yet the real gro wth p ro cess is typically nonlinear, and the surface is non-Gaussian. Indeed, confirmation 



of the scaling given by Eqs. 2^ or 2/I_ in experiments can not be interpreted as conclusive evidence for a self-affine 
geometry: that is a property of the whole ensemble and so requires proper scaling of all moments and correlations, 
not just the second moment. 

1. Quadratic roughness measures 



The quantities (2.2) and (2.4) are quadratic measures of roughness, which we shall also call "linear"; untill recently, 
no other kind was in use. 

The height-correlation function, Eq. ( p.2[ ), is the most standard measure in theoretical discussions, in that "rough- 
ness" is defined by the divergence of this function as its argument r approaches infinity. (Non-monotonic behavior 
of this function has also been used |l5[| to measure the characteristic spatial scale of mounds or other patterns in 
non-self-affine surfaces.) 



On the other hand, the Fourier power spectrum (2^) is central in theoretical derivations but rarely used in experi- 
mental analysis (except for Ref. [p^). This seems to be the best quadratic measure, in that it most cleanly separates 
the contributions from fluctuations on different length scales, and it shows the sharpest knee (in a log-log plot) where 
self-afRne scaling is cut off. 
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Another quadratic measure is the total variance of /i(x) in a in a box of size b, as a function of 5 0Jl|: 

{{h{^)~htf)b, (2.8) 

where hi, = (/i(x))b, and {■ ■ ■)b means the spatial average is only taken over a square of side b centered on xq; this 
variance should be averaged over different choices of Xq. 



C. Fractal properties 



Self-affine surfaces are fractals only in a generalized sense, since the horizontal direction rescales differently from 
the vertical direction. On the other hand, the level set of such a surface (defined as its intersection with a horizontal 
plane) is a fractal object jl9|; see Fig. |l|, below. Different planes of intersection give statistically equivalent level sets, 
since the height fluctuations of a rough surface are unbounded. Level sets consist of contour loops which are the 
connected components. We expect these to be fractal as well, with a fractal dimension smaller then the dimension of 
the whole level set, which is simply the union of all contour loops of the same height. Furthermore, contour loops come 
in all sizes limited only by the system size, and an exponent can be defined that characterizes their size distribution. 
Since contour loops are connected clusters their geometrical exponents are analogous to those defined for percolation 
clusters. 

We will show that the scaling of contour loops uniquely specifies the scaling of the associated self-affine rough 
surface; this will be expressed in formulas giving the geometrical exponents in terms of the roughness exponent a. 
It is somewhat surprising that, by doing measurements solely on the level set, information can be obtained about 
the out-of-plane fluctuations of the surface. For experiments that yield only level-set data without the heights, 
(e.g. freeze- fracture electron microscopy [pO|) our contour- loop analysis is the only route to extracting the roughness 
exponent. 



III. NONLINEAR MEASURES 



In the past the analysis of rough surfaces mostly relied on measures which probed the second moment of the heights, 
such as the height-correlation function or the power spectrum (deflned in Eq. ( |2.2| ) and Eq. ( |2.7] ) below) . But that is 
inherently insufficient to distinguish different growth ensembles or even to verify self-affineness. 

Therefore, to more fully characterize rough surfaces we introduce in this section two new types of non-linear 
measures, i. e., measures that are not linearly related to the structure function of the height field. Non-linear measures 
of the first type (Sec. Ill A) are moments of the "scale-dependent curvature", a modification of the standard height 
correlation function which can identify deviations from Gaussianness of the height fluctuations, in particular the skew 
(up/down asymmetry) at various length scales. (We will compare these to existing non-quadratic roughness measures 



in the Discussion part of Sec. VIII 



Measures of the second type were introduced in Ref. |21| ]; they are distributions of three different geometrical 
quantities deflned for contour- loops (or simply "loops") of constant height, which make up the level sets of the height 
function. These measures are associated with geometrical exponents that characterize contour loops on self-afRne 
rough surfaces: the loop correlation exponent, the fractal dimension of a loop, and the length distribution exponent. 



A. Scale-dependent curvature 



The obvious real-space-based nonquadratic generalization of the height-correlation function is 

(P(x + r)-Mx)]^); (3.1) 

however, this is identically zero on an isotropic surface (and whenever r —^ — r is a symmetry). To escape this problem, 
we observe that h{x -I- r) — ft,(x) is a sort of flrst difference at scale r, and replace it by a sort of second difference. 
Namely, we define the "curvature at x on scale 6" as 

M 

Cb(x) = ['^(^ + ^e™) - '^W] (3-2) 

m— 1 
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where the offset directions {6^} are a fixed set of vectors summing to zero. In our numerical implementation of this 
measure, where {x} is a square lattice, we choose four such offsets related by 90° rotations, pointing either along the 
{10} or the {11} type directions. Those two sets of offsets should give equivalent results (for the same &), provided 
the surface is statistically invariant under rotations in the reference plane. We then define curvature moments (C|) 
for integer powers q. 

The first moment of Cf, is manifestly zero; the second moment is linearly related to the height-correlation function: 



M ^ M 



{Cb{^Y) - M ^ D2ibe„,) "2 2^ D2{b{e,n - e„)) . (3.3) 

m— 1 m.n— 1 



(This is shown by inserting Eq. (3.2) and then decoupling each term of the double sum using the identity (/i^ — 
ho){hn - ho) = ^{ih„^ - hof + [K - h^f - {h„, - ft.„)2}.) 

The higher moments of Cb serve to measure the (possible) deviation of the height fluctuations from the Gaussian 
distribution. For example, if the surface has up/down symmetry h <-^- —h (as all Gaussian surfaces do), ([Cf,(x)]'^) 
vanishes. On the other hand, non-equilibrium grown surfaces often have rounded "hilltops" and sharp "valleys" ; that 
tends to make (C^) > 0, a signature of "skew" in the distribution. Similarly, the fourth moment can also be used to 
test whether the surface is Gaussian, since in that case 

([a(x)]^)/([C,(x)]2)2^3. (3.4) 

For a self-affine surface 

([a(x)]«) ~const|x|«" (3.5) 



follows from Eq. (2.1); of course the coefficient might be zero as is the case for odd q, when the height field has 
up/down symmetry. 

Functions such as (C^) (as a function of b) or Dq{r) (as a function of r) can also be used as "spectra" of the height 
fluctuations for non self-affine surfaces. That is, differences in the behavior of the function in different ranges of b or 
r reveal qualitative differences of the surface morphology on the corresponding length scales. 

In principle, q-th or der moments may scale with well-defined exponents aq, yet the surface is not self-affine since 
aq qa violating Eq. ( |3.5| ); this is called a "multifractal" or, more precisely, "multiaffine" surface Slow transients 
of multiaffine behavior (up to ~ 10^ steps ind = 1 + 1 and ~ 10'^ steps ind = 2+1) have been seen recently in numerical 
simulations of growth models |23| (which, however, are believed to be asymptotically self-affine). The analogous higher 
order structure functions are a central issue in turbulence, where the violation of self-affine (Kolmogorov) scaling is 
well established and is associated with intermittency of the velocity field fluctuations . 

The scale dependent curvature can be contrasted with Krug's height-difference moments p^ , 

Dq{r)^{\h{^ + r)-h{^W), (3.6) 
a na tura l generalization of the height-correlation funct ion using an absolute value to avoid the trivial cancellation in 



Eq. (3^). Das Sarma and collaborators [p5[|23| ] used (3^) to test for multi-affine behavior (whereby the 1/q pow er 



of the q moment scales with exponent depending on q, unlike the simpler self-affine case). For odd q, Eq. ( |3.6D is 
insensitive to the up-down symmetry (or lack thereof) since it is nonzero anyhow. Our "curvature" seems to be the 
simplest function that detects the skew locally. 



B. Fractal dimension of contour loops 



For the remainder of this section, we must define the loop ensemble. Consider a contour plot of a rough surface 
with a fixed spacing A between heights of successive level sets. We take it to be an arbitrary constant much smaller 
than the typical (r.m.s.) fluctuation of ft.(x). The value of A does not affect our exponents and we need to consider it 
explicitly only in the arguments of Sec. IV A; in other places we may implicitly scale /i(x) such that A = 1. In STM 
images of rough metal surfaces A is usually the height of a single step on the surface. 

The contour plot consists of closed nonintersecting lines in the plane that connect points of equal height, which we 
call contour loops (see Fig.^ below). Every random-surface conflguration maps to a configuration of contour- loops; 
when the probability weights of the respective configurations are taken into account, this deflnes a mapping of the 
random-surface ensemble, to the contour-loop ensemble. The contour loop ensemble arising from self-affine random 
surfaces is (we shall argue) self-similar; the loops are connected clusters that can be studied using scaling, just as 
(critical) percolation clusters have been analyzed in previous work Q . 
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For every contour loop in the loop ensemble we define a loop length s and a loop radius R. In all the examples we 
study the heights are defined on an L x L square lattice with lattice constant a. The loop length is measured with a 
ruler of length a while the loop radius (really a diameter) is defined as the side of the smallest box that completely 
covers the loop; see Fig. |l8|. 

In the loop ensemble we define a joint distribution n(s, R) (independent of the contour spacing A) such that the 
number of loops with length in (s, s + ds) and radius in {R, R + dR), per unit area, is 

A-^il{s,R)ds dR. (3.7) 

The factor A^^ has the obvious significance that if one halves the contour spacing, one has twice as many contours. 
Assuming that the loop ensemble is scale invariant, we expect that n(s, R) has a scaling form 

n{s,R)^ s-yUis/R^'f) . (3.8) 

Here Df is the fractal dimension, and y is simply related to the length distribution exponent r, which we define in 
the next section. 

In practice the exponent Df is measured by the scaling relation 

{s){R)^R^f, (3.9) 

where 

pR+SR. j-R+SR 

{s){R)= sn{s,R)dsdR/ n{s,R)dsdR (3.10) 

Jr Jr. 

is the average loop length for loops whose radius falls in th e in terval (i?, R + SR), 6R <C R- The scaling in Eq. (|3.9|) 



follows immediately from the assumed scaling form in Eq. (3.i 

The dimension defined in Eq. ( |3.9D is really the scaling dimension of the loop length, i.e., it defines the relation 
between bigger and smaller loops in the distribution. On the other hand, the proper fractal dimension (either the 
Hausdorff dimension Dh or the self-similarity dimension) refers to the relation between bigger and smaller pieces of 
the same loop. Thus Dh is defined by s a~^" , i.e. how the loop length scales with the ruler size. When the 
contour- loop distribution is self similar (as we shall assume), the two kinds of dimensions are equivalent. 



C. Loop length distribution exponent 



We define the loop number density P{s) so that A^^ P{s)ds is the total number of loops, per unit area (measured 
in sites), with lengths in (s, s + ds); a related distribution of loop lengths, P{s), is defined such that A~^P{s)ds is the 
number of loops passing through a fixed point (say the ori gin) with lengths in the range (s, s + ds). In lattice models 
(including our numerical examples in Sections |V|, VI, and VII), s is an integer and P{s) is essentially the probability 
that the loop has length s. 



From comparison to Eq. (3.7) it is obvious that 



P{s) 



n{s,R)dR. 



Since the total number of sites along a loop is equal to its length s we have 

P(s) = sP(s); 

the additional factor of s is because each site could be the origin in the definition of P{s). 

Assuming that the loop ensemble is scale invariant we can define the length distribution exponent r by 

P(s) - s-(^-i); Pis)r^s~\ 



(3.11) 



(3.12) 



(3.13) 



This is t o hold for l arge contour loo ps, i. e. those of radius much bigger then the microscopic scale a. Indeed, inserting 
Eq. (U) into Eq. (|3.1l| ) gives Eq. ( |3.13| ), with 



y - T + 1/Df 



(3.14) 
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On the other hand, we could also define n(i?) such that A ^n(R)dR is the total number of loops, per unit area, 
whose radius is in the range {R, R + dR) . Obviously 



n{R) 



n(s, R)ds. 



(3.15) 



Doing the integral and then eliminating y using (3.14) gives 

n(i?) - l/i?i+^^("-i) (3.16) 
We would have obtained the same result more quickly (and more dubiously) had we assumed a strict relationship 



between radius and length, s — (const) R^ f , rather than write (3. 



D. Loop correlation function 

The loop correlation function G(r) measures the probability that two points separated by r lie on the same contour 
loop. This correlation function is non-local, for the connectedness of the two points depends on every site on the 
portion of loop between them. This loop correlation function should be distinguished from the level-set correlation 
function which simply measures the probability that two points separated by r are at the same height. For the loop 
correlation function to be well defined the contour lines are considered to be of finite width given by the microscopic 
scale a. Due to rotational symmetry of the loop ensemble, G(r) depends on r = |r| only, and for large separations 
(r ^ a) we expect it to fall off as a power law: 

G{r) ^ -L , (3.17) 

This equation defines the loop correlation exponent xi which is at the heart of the scaling theory of contour loops 
developed below. 



IV. SCALING RELATIONS 



In this section we derive scaling relations among the roughness exponent a, and the three geometrical exponents 
- Df, r, and xi - associated with contour-loops and defined in Sec. III. These formulas are corollaries of the self- 
affineness of the rough surface, Eq. (2.1). Furthermore, for growth on an initially flat substrate the heights will be 
uncorrelated beyond a certain time-dependent length scale and the large contour loops are best modeled as hulls of 
percolation clusters. This implies a crossover to a different set of expon ents as w orked o ut in Sec. IV E. The scaling 
relations - including the finite-size and finite-time forms in Sec. IV Dl and Sec. IVEl - will serve as a useful tool 



for analyzing the surface morphologies obtained from numerical simulations and in experiments (see sections VI and 
Vl| . 

There are three stages of the main derivation. First, we establish a relationship between the self-affine exponent a 
and the loop-size distribution exponent r; it is analogous to the hyperscaling relation among percolation exponents. 
Second, we fi nd a sum rule (analogous to the susceptibility sum rule) relating the loop correlation exponent xi of 
Section HID and the loop-size distribution exponent r. Third, we present a conjecture that the loop correlation 
exponent has a value xi = 1/2, which is super-universal in the sense that it is independent of a. (This conjecture is 
supported by an exact calculation of xi in the extreme cases, i.e. a = (equilibrium rough case) and a = 1.) Finally, 
these relations taken together yield formulas for Df and t (Eq. (4.13)) as a function of a. 



A. Hyperscaling relation 



If we parametrize a loop as l(s), where s is the arc length as measured by a ruler of length a, then after the rescaling 
given by Eq. (2.1) it is mapped to, 



l(s) -> b-'lib'^fs) 



(4.1) 



This scali ng proper ty of the contour ensemble justifies the power law dependence of G(r) on r and P{s) on s, in 
Eq. (|3.17D and Eq. (|3.13|) respectively. 
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In writing Eq. (4.1), we made a nontrivial hypothesis that the contours of the height function obtained by coarse- 
graining a given reaUzation of /i(x) are statistically the same as the coarse-grained version of the contours of h{x). We 
know of no coarse-graining procedure for the height function which assures that the contours will stay the same. It 
will happen that, near a saddle-point of ft.(x), two loops (both of height /iicv) approach closely, but the coarse-graining 
shifts the height of the saddle-point across ft-iov so that the coarse-grained versions of the loops coalesce into one loop. 
Whether this phenomenon makes a relevant contribution to our scaling relations depends on the frequency of close 
approaches pTj . 

To determine the scaling of n(i?) first apply the rescaling Eq. ( pj| ) to each configuration of h{r); this maps the 
contour ensemble to a new contour ensemble with rescaled contour interval A' = 6^" A. The total nu mber of contours 
with radii in the range {R, R + dR), in a box of side L, is L'^ A~^n{R) dR, by our definition in Sec. Ill C . Since the 



contours are mapped 1-to-l (according to the hypothesis in Eq. (4.1)), we can equate this with the number of new 
contours in a box of side L/b and of radius in {R/b, R/b + dR/b), which is {L/b)^A'~^n'{R/b)dR'/b. On the other 
hand, by self-affineness the new height ensemble, is statistically identical to the original one (for large R); this holds 
as well for the new contour ensemble, thus fi'{R) = n{R). So we obtain n{R/b) — b^^°'h{R), which implies the scaling 
behavior 



fi{R) - R 



-3+a 



Equating Eq. (O) and Eq. ( 3.16 ) leads to the first scaling relation (called "hyperscaling" ) 

DfiT-l) = 2^a. 



(4.2) 



(4.3) 



This scaling relation has been derived previously by Huber et al. p8| in a slightly different context, and in a somewhat 
different form by Isichenko and Kalda |^ . Unlike the usual hyperscaling relation for percolation clusters which can 
be derived from the assumption that the number of large clusters does not grow with scale of observation, here that 
number grows as a power with exponent a pffl. 



B. Sum rule 



A second scaling relation can be derived from a sum rule. To start off, let's separately consider the loop correlation 
function for different loop sizes s. Let Gs{v) be the probability that point x -I- r is on the same loop as x, given that 
the loop has length s. In light of the self-similarity of the loop ensemble, it is reasonable to assume 



Gs{v)^s"^\v\-^fGs{rls^/^f) 



(4.4) 



where m and a are as-yet undetermined exponents, and /gs() is a scaling function. The reason we must scale r by 



is that this is the typical diameter R of the loop (s ~ R ^ ). 



Now, the sum of G's(r) over all lattice points is the expectation of the total number of points in the loop, which 
was given to be s, hence (substituting from (L4)) 



d'^rGsir) ~ s 



{2-a)/Df+m 



which gives one relation between the exponents a and m introduced in (4.4): 

2 - a = Df{l - m) 



(4.5) 



(4.6) 



On the other hand, the total loop correlation is the integral of Gs(r) over the loop distribution fmiction P{s) given 
by Eq. (p^, thus 



G(r) = / dsP{s)Gs{r) - r-''(r^O 



-a( Df\m,+2-T ^ Dy(3-T)-2 



(4.7) 



where ( |4.6| ) was used to eliminate both a and m in the result. Equating the exponent of G(r) in (4/7) to the one 
defined by (3.17), we obtain the scaling relation 



Df{3-T) ^2-2x1 



(4.8) 



The above scaling relations, Eq. (4.8) and Eq. ( |4.3| ), can be combined into expressions (which were originally 
presented in Ref. |^l|) for the fractal dimension Df and the exponent r: 
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Df = 2-xi- a/2 



(4.9) 



r-l = ^ " , . (4.10) 

2 -XI- a/2 ^ ' 

The first scaling relation is reminiscent of the relation 

D^2~a (4.11) 



due to Mandelbrot ||T^. The important difference is that Eq. ( 4.11 ) gives the fractal dimension D of the level set of a 



random seif-affinc surface, and not the fractal dimension of a single contour loop. (We emphasize this point because 
there has been some confusion in the literature where the two dimensions have been equated.) 

Olami and Zeitak considered the same loop ensemble, but mostly focused their attention on the "islands" 
contained in the loops rather than the contours; their "t" exponent (which we call tzq) refers to the distribution of 
island sizes. They derived a formula tzq = 2 — q;/2 (in our notation). It is easy to show 2(rzo — 1) = Djij — 1) - the 
"2" here is the fractal dimension of these islands [32|; upon inserting this conversion, their formula turns out to say 
Df{T — 1) = 2 — a, which is the same as our Eq. ( |4.10| ). 



C. Loop correlation exponent 

Now we turn our attention to the contour correlation exponent, and we conjecture that 

XI - 1/2 (4.12) 

is super-universal in that it is independent of a. 

In the case of an a = Gaussian surface, we know xi = 1/2 exactly for a solvable statistical-mechanics model of 
contour loops, equivalent to the critical 0{2) loop model on the honeycomb lattice Details are in Appendix 
By invoking universality this is valid for all logarithmically rough random Gaussian surfaces. 

The exact value of xi can also be determined for a — 1. Namely, the fractal dimension (D f) of a contour loop must 
satisfy Df < D since it is a subset of the level set, which has dimension D = 2 — a = 1; Eq. ([4.11 ). On the other hand 
Df > 1 since a loop has topological dimension one. From these inequalities we conc lude that for a = 1 the fractal 
dimension of a contour loop is Df = 1. This in turn leads to xi — 1/2, from Eq. 



The validity of conjecture (4.12) for general a has been checked, to date, only through the numerical simulations 
reported in Section M and in numerical simulations of Zeng et al. [ p4[ . 

Since xi — 1/2 for a — and a = 1, a proof of monotonicity of xi with a would suffice to establish the conjecture. 
Even that is very difficult owing to the non-local definition of the loop correlation function. 

D. Combined scaling relations 

Equipped with the (super-universal) conjectured value of the loop exponent xi — 1/2, and the scaling relations. 



Eqs. (4.9) and (4.1C), we find the following formulas for the geometrical exponents of contour loops of a self-affine 

3 — a 



surface with roughness exponent a: 



r-l^i^^ (413) 

3 — a 

These relations form the basis of the contour loop analysis of rough surfaces, which we implement in the following 
sections. 

Our formula for D f differs from the one proposed by Isichcnko | p5[ | 



£,Isichcnko ^ ]^ ^ ^ (4.14) 



which was derived from an approximate "multiscale" analysis. We note that the formula for I?/ in Eq. (4.14) gives 
the wrong result in the a = case, where Df — 3/2 is exact | |3^ ]. 
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1. Finite-size scaling 



For realistic rough surfaces the self-afRne scaHng wiU be cut off at large lengths either by the correlation length or 
the system size. 

In the case that self-afhne scaling is cut off only by the system size L, we can extend the power laws derived above 
for the average loop length, the size distribution of loops, and the loop correlation function, into scaling forms: 

= i?^//,(i?/L) 
P(s,L) = s-(--i)/p(s/L^/) 

G{v.L)^\v\-^-^^f{\v\/L) . (4.15) 

In the case that the self affine scaling is cut-off by a finite correlation length ^(t) < L, our three contour-loop measures 
will display crossover effects to a different set of power laws; we turn to this problem next. 



E. Percolation crossover 



For a surface roughened by growth self affine scaling is expected to hold only up to a finite correlation length ^{t) 
growing with time as Eq. ( ^.3[ ). At early enough stages of growth (i.e. while ^(t) < L) the statistics at scales beyond 
^(i) depend on the initial state. Then the contour loops of the surface will also exhibit crossover behavior where loops 
whose linear size (as measured by the radius R) is less than the correlation length will scale according to the formulas 
derived above, while the large loops will exhibit scaling with percolation exponents. 

Say the initial surface is flat (which we assume henceforth) . Then it turns out (see Appendix ^ that the contour 
loops at scale R> £^{t) are boundaries of percolation clusters. Although this new contour loop ensemble corresponds 
to a noK-self-affine surface, it still exhibits scaling and we derive its three loop exponents in terms of known percolation 
exponents. In some cases, it turns out, the exponent values from the percolation regime and from the self-afRne surface 
are not so different; thus a careless analysis might yield spurious exponents. 

It is easy to see that at r > ^, we can model the actual heights (not height differences) as statistically independent, 
since (by definition of the correlation length) the distance ^(t) is the farthest that an event can influence another in 
time t. 

In appendix we derive the geometrical exponents 

Df^p = 7/4 = 1.75, Tp = 18/7 = 2.571, 2xup = 5/4 = 1.25 (4.16) 

which apply to loops at scales larger than ^{t) (the percolation-regime scaling). These exponents are the same for any 
a. An important corollary is that power-law scaling in the loop analysis is not necessarily a signature of self-afRne 
behavior. Indeed, most real surfaces never reach a clear self-affine regime, hence their loops are probably in the 
percolation regime. 



1. Finite-time crossover scaling forms 



Th e complete crossover between the self-affine and percolation regimes is described by scaling forms parallel to 
(4.15). First consider the height structure factor. Since, as noted above, the heights are independent, their (spatial) 
power spectrum is flat in Fourier space: S'(q) ~ const, for |q| < l/^{t)\ on the other hand, for |q| > 1/C(0 the surface 



has already developed a self-affine state so Eq. (2.7) does hold. The two behaviors should be combined via a scaling 
function /s(): 



5(q;t) = |q| 



-2(l+a) 



/s(qt'/^) 



(4.17) 



see Fig. 1(a). 

Thus, at times t such that ^(t) ^ L, the behaviors (3.9), ( 3.13| ), and (3.17) are generalized to 



(4.18) 
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In each case, the scaling function is unity for argument zero, while for argument large it s cales as a power law needed 



to give the correct exponent for the percolation regime, as calculated in Appendix B2. The t dependence for the 



prefactor of each percolation-regime power law is given by the requirement to patch the above two dependences 
together when the scaling-function argument is of order unity. 

Figures^ (b)~(d) illustrate the shapes of the three loop measures. Notice that the "knee" around r = ^{t) appears 
more strikingly in the Fourier analysis than in any of the loop analyses. Although the percolation- regime and self- 
afhne exponents have fairly similar values, the difference grows larger as a gets larger. The crossover is evident in our 
simulated Gaussian data (see Fig. [lo|) . 



V. SIMULATION: GAUSSIAN RANDOM SURFACES 



Here we test the validity of our scaling relations and the effectiveness of determining a from contour loops, under 
the controlled circumstances provided by computer generated surfaces with known a. The surfaces we construct are 
self-affine with Gaussian fluctuations of the height. 



A. Construction 



Random Gaussian surfaces are generated numerically as an L/ax L/a matrix /i(x) of real- valued heights associated 
with the vertices {x} of a square lattice of size L, with lattice constant a. A particular realization of /i(x) is given by 
Fourier transforming /i(q) where the wave- vectors q take their values in the first Brillouin zone [— tt/o, tt/o] x [tt/o, tt/o] . 
Each Fourier component h{q) is an independent Gaussian random variable with a q-dependent variance given by 

(IMq)P> = . (5.1) 

For > a > 1 surfaces generated in this way are self-affine and rough, with a roughness exponent a. 

The a = case of random Gaussian surfaces is familiar as: (i) the equilibrium-rough surface (compare Eq. ( |A1[ )), 
(ii) the surface in the Edwards- Wilkinson model, and (iii) the Coulomb gas representation of a two-dimensional critical 
model (see appendix^. The case a — I appears in the MuUins-Herring (diffusive relaxation) model (a = 1) of 
non-equilibrium surface growth |]3^ . 



1. Comparison to other studies 



A popular algorithm for generating self-affine surfaces is "random midpoint displacement, with random successive 
addition" ||3^,|l^. This method iterates a step in which, starting with a self-affine surface on a coarse grid of lattice 
constant 2a, one generates heights on a new grid of lattice constant a by interpolation, and then adds to them random 
increments proportional to a". Such an ensemble need not be Gaussian or have up-down symmetry, but commonly 
does [0. We note that (i) the variance of a site's height (relative to the initial flat surface) depends on what iteration 
that site appeared, i.e. on how many times 2 can be divided into the site coordinates; (ii) height-difference correlations 
do not always grow with distance (they are smaller between two sites that appeared in early iterations) and they 
have the anisotropy of the lattice even at large distances. We believe our Fourier construction of self-affine surfaces 
(Sec. VA) is preferable because the resulting ensemble is (i) spatially homogeneous and (ii) isotropic, on scales beyond 
a couple of lattice constants. 



B. Curvature measurements 



We measured moments (C™) (m = 2,3,4) of the scale-dependent curvature, as defined in Sec. Ill A, for Gaussian 
surfaces generated by the Fourier method described above. This data (in Fig. |^) is a kind of check on the Fourier 
method since the mean over an infinite number of samples can be computed analytically. 

Self-affine scaling is evident on the log-log plot of the even moments in Fig. || (upper plot). The roughness exponent 
a is obtained as the slope of a straight-line fit to the plot, as shown in Table Ideally, the slopes of the and 
log- log plots should be 2a and 4a with exactly the input a values used in constructing the random surfaces. This is 
spoiled somewhat in practice by discrete-lattice effects for 6 < 3 and by finite-size effects when b > L/A. Furthermore, 
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(C'b ) / (Cfc )^ should be exactly 3 for every b value, even those for which the power-law dependence on b fails, since this 
is true for any Gaussian random variable. Indeed, the measured ratio is close to 3. 

The third moment of Cf, is shown in Fig. ^. Independent of a, (C^) is roughly zero, as expected for a random 
Gaussian surface which posses a /i — > —h symmetry (i.e., the valley bottoms and the hill tops are equivalent for a 
Gaussian surface). 



C. Loop measurements 



The primary motivation for our Gaussian surface simulations was an initial test of the scaling predictions for the 
contour loop exponents from Sec. IV. A contour plot of a sample surface configuration for a = 0.4 is shown in Fig. 
(A similar plot for a — was published in Ref. |21|.) 



1. Measurement procedure 



In a single run, which would typically take 10 minutes on an Sun SparcS workstation, 25 surfaces of specified 
roughness a were generated. For each surface typically 400 points were chosen at random, and through each point 
a contour loop was constructed using the loop finding algorithm as explained in Appendix |c[ While each loop was 
being traced points along the loop were used to evaluate the loop correlation function G{r). For each contour loop 
its radius and length were measured and used to determine the length distribution of contour loops (-P(s)) and the 
average loop length (s) as a function of the loop radius R. 



2. Results 



In order to measure the geometrical exponents Df, t, and xi we plotted the data for system size L = 512 on a 
log-log graph and performed least-squares linear fits. Data was selected for fitting from the range in which a well 
developed power law was observed; see Figures ^, ||, and ^. The results are given in Table |. We find excellent 
agreement between the predictions of the scaling theory and the measured geometrical exponents. In particular, note 
that the simulations confirm the super-universal nature of the loop correlation exponent xi = 1/2. 

The loop correlation function G{r) has a a size dependence which biases a direct fit to the exponent 2xi; finite- 
size scaling (see below) partially overcomes this systematic error. Our theory (Sec. IV C) indicates that G(r) has 
a universal exponent 2xi — 1; in fact, as shown in Fig. ^, G{r) itself appears practically independent of a. Closer 
examination reveals that the coefficient in G(r) ~ 1/r decreases slightly as a grows. Furthermore, the fitted values 
of 2xi (see Table |) decrease a bit with a, which we attribute to the systematic error just mentioned, combined with 
the small a-dependence of the shape of the "knee" in the finite-size behavior of Fig . ^ There is no indication in 
the extracted 2xi values of any non-monotonic dependence on a; as shown in section [VC, monotonicity of xi(a) is 
sufficient to prove 2xi — 1, independent of a. 

A better measure of th e geo metrical exponents was obtained from a finite size scaling analysis of the data. Using 
the scaling forms in Eq. ( 4.15| ) we produced data collapses ("scaling plots"). Sample data for the a = 0.4 case are 
given in Figs. |^(a) and ^(a); the data collapse is shown in Figs, ^(b) and ^(b). From the loop-size distribution plots 
like Fig. 0(b), we extracted both the exponents Df and r — 2. Similarly, we obtained 2xi from the loop correlation 
function Fig. ^(b); in this case we don't fit another exponent since r obviously scales as L^. (We did not carry out 
finite-size scaling of the (s) versus R plots such as Fig. |^, since there was no obvious change in the slope as a function 
of R/L.) The geometrical exponents giving the best data collapses are reported in the "FSS" columns of Table |. 
The reported uncertainties were estimated by the interval over which changes in the exponent value did not visibly 
worsen the data collapse. 

Note in Table | how the finite-size scaling exponents agree better with the scaling theory of Sec. ^ than the 
exponents obtained from "direct" fitting of the data to power laws. The discrepancy becomes more obvious at larger 
values of a. We infer from this that Gaussian surfaces with a large value of the roughness have more pronounced 
finite size effects which lead to an overestimate of Df and r. This is of relevance to experimental data where the 
system size is typically not a tunable parameter, and the geometrical exponents are necessarily measured using the 
direct-fit method. 
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3. Relation to a previous simulation 



Numerical measurements of the fractal dimension of contour loops have been done by Avellaneda et a l. [p9| . They 
found Df — 1.28±0.015 for an a = 0.5 surface, which is close to the predicted value Df — 1.25, from Eq. (^.13| ). They 
also measured the combination Df{T — 2) (their "a") which describes the scaling of the probability that a loop passing 
through a fixe d point has a radius larger than p, with p. (We evaluate this quantity by integrating n(s, R) = sn{s, R), 
from Eq. (|3.8| ), over all s and for R > p). The numerical result they quote, Df{T — 2) = 0.21 ±0.017, is in fair 
agreement with our prediction -D/(t — 2) = (1 — q;)/2 = 0.25 for a = 0.5, which follows from Eq. ( [1.13[ ). 

Wagner et al jssj simulated a form of invasion percolation where the threshold pressures have the form of a self-afRne 
surface. Hence the perimeters of the invaded clusters are the same as the contour lines of the surface. They claim that 
the perimeter dimension is "consistent" with Isichenko's formula, our Eq. ( 4.14 ), but do not quote an error; perhaps 
their precision was suc h th at the prediction of (4.14) could not have been distinguished numerically from the one we 
beheve to be correct, ( [4.13| ). 

Ref. | |3^ ] also mention measuring a behavior r ^ with 7 w 0.9 for the correlation between successive filled sites. 
If this were simply a correlation of two randomly chosen points along a perimeter, it would be identical to our loop 
correlation functions, G(r) or Ggir) defined in (3.17) or (4.4); in fact, the filling process would appear to depend on 
correlations of the surface gradient and might have a somewhat different exponent. 



D. Surfaces with a finite correlation length 



To test the percolation analysis of self-afhne rough surfaces with a cutoff, as derived in Appendix ^ and summarized 
in Sec. IV E| , we performed curvature and loop measurements on Gaussian surfaces with a correlation length ^. The 



correlation length is incorporated in the Fourier method of generating Gaussian surfaces by changing the variance of 



/i(q) in Eq. (5.1) to: 



The effects of the cutoff are summarized in Fig. |l^, which should be compared to the theoretical prediction of 
Fig. 1^. The curvature and loop data shown in the figure are for system size L = 512. 

The second moment of the curvature displays self-affine scaling with roughness a = 0.4 up to a length scale set by 
^g, and beyond this scale it levels off; see Fig. p^a). We checked that the third moment of the curvature vanishes, as 
expected since the height fluctuations are still Gaussian, while the fourth moment follows affine scaling up to roughly 
the same correlation length as the second moment. 

The loop measures exhibit distinct crossover behavior, as seen in figures Fig. |l^ b) through d). For loops whose 
radius is smaller than the correlation length, which is here ^ w 20, we find values of the geometrical exponents 
consistent with those extracted previously for a = 0.4 random Gaussian surfaces. For loops whose linear size exceeds 
the cutoff, scaling consistent with the percolation analysis is found. The actual numerical values extracted by fitting 
the 20 < i? < 200 data to a power law are somewhat larger than expected (2x/ — 1.46(6), Df — 1.7(1), and 
T = 2.63(1)) which we attribute to finite size and/or crossover effects. To check this we also simulated a Gaussian 
surface with completely uncorrelated heights, i.e., with ^ = 1 and system size L = 512, for which we find (by the 
direct- fit method): 

2x1 = 1.26(3) , Df ^ 1.70(2) , t = 2.565(10) (5.3) 



in good agreement with Eq. (4. IE) 



VI. SIMULATION: NON-EQUILIBRIUM GROWTH MODEL 

In this section, the linear and nonlinear roughness measures of Sec. [II, which in Sec. ^ were tested on artificial 
Gaussian random surfaces, are now applied to growth-roughened surfaces produced by a simple random deposition 
model, the well-known "single-step model" . Our results are in support of the view that the single-step model produces 
self-affine morphologies. 
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A. The single-step model 



We implemented the "single-step model" pO|-|4^ in d = 2 + 1 dimensions . (More details on this model are found 
in Sec. Ill F of Ref. |^^, or Sec. II A of Ref. ||42||.) There is one control parameter The allowed configurations 
are just those of the BCSOS model: each site of a square lattice has an integer-valued height and neighboring heights 
must differ by ±1. The Monte Carlo rule is that in each time step a deposition event occurs with probability p+ or 
an evaporation event (inverse of a deposition event) occurs with probability 1 — once it is decided which type of 
event occurs, a site is picked at random among those sites at which that event is allowed Q]. 

We begin by an overview of the theoretical expectations. Up-down symmetry switches p+ ^ 1 — thus we need 
only report data for < p+ < 0.5. The case p+ = 0.5 is special as the dynamics satisfies detailed balance. This should 
produce an equilibrium-rough interface, namely the BCSOS model with all configurations weighted equally p5| . This 
interface, at long wavelengths, is described by the Gaussian model of Sec. ^with a — (Edwards- Wilkinson behavior). 

On the other hand, the growth model for p+ 7^ 1/2 is believed to asymptotically belong to the Kardar-Parisi-Zhang 
(KPZ) universality class It has been proposed that a = 0.4 exactly for the 2+1 dimensional KPZ model 

p6| ; however, finite-size effects, small simulations and naive fits systematically underestimate it as a ~ 0.38 
The KPZ behavior should be clearcut when p+ is close to 1, but otherwise a crossover from initially Gaussian to 
asymptotic KPZ behavior is expected, which will be slow (as a function of time or system size) if p+ is close to 1/2. 

It turns out, in our numerical results (below), that p+ = 0.5 indeed shows Gaussian behavior and p+ = 0.1 shows 
KPZ-like behavior, but p+ = 0.3 consistently resembles p+ = 0.5, at the sizes we could simulate (i.e. up to L = 128). 
We attribute this to above-mentioned crossover from initial Gaussian behavior. 



B. Simulations 



Starting from a flat surface, we ran the simulation (for systems of 128 x 128 sites) for 2000 Monte Carlo steps (MCS) 
per site to equilibrate and then took data for a period of 1200 MCS/site; one such run took 10-15 hours of cpu time 
on a RISC-6000 workstation. The standard-length runs (for size L — 128) appeared to be insufficiently equilibrated 
for p+ = 0.1, since they failed to collapse on finite-size-scaling plots with smaller systems. Therefore we performed 
one run for L — 128, p+ = 0.1 with 12000 MCS/site equilibration and 10000 MCS/site for data collection; this is the 
run reported in our results. In all other cases, we believe the run time was adequate, since much shorter runs showed 
no gross differences. We performed about four runs for each value p+ = 0.1,0.3, and 0.5, verifying the symmetry 
p+ ^ 1 —p+. (All measures are the same, apart from a change in sign of (C^).) Only one of the ~ 4 runs was selected 
to be fitted and plotted here; the data sets presented as e.g. p+ — O.l are actually p+ — 0.9 in some cases. 

Once every 100 MCS/site (during the data-collecting portion of a run), we performed a measurement step on the 
surface. The Fourier transform was taken of h{r) using a fast-Fourier-transform routine, but (|/i(q)p) was accumulated 
only for q values along the (1,0), (0,1), (1,1), and (1,-1) directio ns. A lso, in each measurement step 100 contour loops 
were traced out from random initial points, as described in Sec. VC and Appendix |c[ Statistics were accumulated of 
the loop's radius R and its number of sites (length) s, but not the loop correlation function. 



C. Fourier and curvature results 



The single-step model is the only one simulated in this paper for which we evaluated Fourier spectra, which are 
plotted in Fig. |l^; the log-log plot should have a slope —2(1 -I- a) so a can be extracted from a linear fit (as in Table ||). 
Notice how the spectra are completely isotropic with respect to the lattice directions. 

The scale-dependent curvature moments were not evaluated during the runs, but were computed only from the final 
surface from each run (henc e thei r statistics are much worse than for other measures reported here). We computed 



(C6(x)™) as defined in Sec. [II A, for m = 2,3,4, as a function of the offset b in the definition of Cb as a discrete 



Laplacian. The results for m = 2 and 3 are plotted in Fig. ^ (All figures of the SSM are from the largest system 
size, L — 128.) Fig. |l^(a) does not show well-defined power laws. The curve for p+ = 0.5 shows a smallish apparent 
slope, 2a = 0.3(1), and a downwards curvature which is plausibly consis tent with the expected logarithmic behavior, 
just like that of the usual height-difference function D2{h) (recall Eq. (3.3). The (C^) curve for p+ = 0.1 shows a 
larger slope, 2a ~ 0.6(1), consistent with KPZ scaling. The a values in Table were extracted from fits to (C^) 
plots. Slopes from plots of the (C^) moments (not shown) are consistent with 4a for the a values in the table. 
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D. Loop analysis 



We analyzed the loop ensemble to plot the mean loop size as a function of its radius (Fig. |l^ and the cumulative 
loop-size distribution P>(s) (Fig. p^ ). (Note that s must have even values' thus it is necessary to divide the nonzero 
values by 2 in order to properly estimate P{s), which was assumed (in Sec. |lj and Sec. [V) to be a smooth monotonic 
function.) 

As with the Gaussian simulation of Sec. ^ the slopes of straight-line fits to log-log plots of these data, give estimates 
of Df and r — 2; they are tabulated as "direct" in Table Alternatively, we used data like Fig. 1^ from smaller 
sizes L = 32 and i = 64 (as well as L = 128) to produce scaling plots analogous to Fig. [t] (these plots not shown), 
extracting the "FSS" data in Table |^ (We do not report on 2xi since we did not evaluate the loop correlation G{r) 
in our SSM simulations.) 

It is interesting to compare the four different measures of a included in Table ||. Those from (|/i(qp) seem to have 
the smallest statistical errors (and the most sensible values). The closely related (C^) result is expected to be worse, 
not only because the statistics are poor in our implementation of the SSM simulation (see above), but also because it 
uses /i(x) values from more widely spaced x and is therefore more sensitive to the system size. 

The next best method seems to be the Df loop analysis from ((s),i?) plots; curiously, it appears that the Dj fits 
show smaller run-to-run fluctuations than the (|/i(qp) fits. As also observed in the Gaussian runs (Sec. 0), the P>(s) 
analysis showed more obvious finite-size effects; direct fits to t are unreliable and only finite-size scaling plots give 
reasonable results. 



VII. ANALYSIS OF EXPERIMENTAL DATA 



In this section we test our nonlinear measures against experimental scanning tunneling microscopy (STM) data sets. 
Rough metal surfaces grown under several conditions are believed to devel op a m orphology with self-afRne scaling, 
but only up to a time-dependent correlation length ^{t) as discussed in Sec. IV E . 

The most detailed analysis was done for the vapor deposited Ag surface on a quartz substrate of Palasantzas and 
Krim |^^. We obtained a 400 x 400 height array corresponding to an STM image of a 702nm thick Ag surface and 
performed curvature and loop measurements. Note that all the results quoted below are from a single height profile. 
All in-plane lengths will be measured in units of the grid of this data, which is 1.625 nm. 

We also report briefly (Subsec. VII C) a less thorough analysis of an STM data set from a different, but still 
self-affine, growth regime showing KPZ scaling. 



A. Quadratic measures and curvature moments 



Palasantzas and Krim p9| | originally evaluated the roughness exponent 

a = 0.82(5) (7.1) 

from a fit to the standard (quadratic) height correlation function D2{r) defined in Eq. (|2_ 

That correlation is similar to our second curvature moment (C^). This quantity shows a power law dependence on 
the scale &, up to a correlation length which was estimated to be ^ = 25(5); see Fig. ttq(a). A linear least squares fit 



of the data with 6 < ^ gave 2a ~ l-^Q), agreeing (as expected) with ([7.l|). 

The third moment (C(,(x)^), Fig. [I5|(b), shows distinct non-Gaussian behavior, as expected for non-equilibrium 
growth. It reaches a maximum at length scale 6 « 23 which correlates well with ^. This indicates a morphology 
consisting of grains of typical size ^ that are rounded at the top. Such a morphology is clearly seen in three-dimensional 
renderings of the STM data in Ref. |^ , or the gray-scale image in Ref . |5^ . 



1. Kleban's nonlinear measure 



Recently, Kleban et al defined a non-quadratic roughness measure rather different from any of those mentioned in 
Sec. m. First, for every scale b they constructed a smoothed version i?f,(x) of the height function, as the average of 
/i(x') for x' in a 6 X 6 square centered at x (alternatively by convolving with a Gaussian weight function of width 
b.) Then they calculated the histogram P{Hb) of i?f,(x) values for x ranging over the entire sample, and the skew 
moment of this distribution. (Of course, P{Hh) is defined for a rough surface only when a finite-size or finite-time 
cutoff is present.) 
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When applied to the STM data of Palasantzas and Krim |^^, the distribution of Hi,{x) appears Gaussian when 
b = 0, i.e. for the raw data. That is rather mysterious, since the surface certainly lacks up/down symmetry: it consists 
of deep, narrow crevasses and rounded hills. Each crevasse should contribute to a long tail on only the h < h side 
of the distribution, while each hill contributes a peak and a sudden drop to zero on the h > h side. However, the 
fluctuation in height from hilltop to hilltop smears out that sharp feature, making a spuriously symmetric distribution. 
As pointed out in Ref. their smoothing of h{x.) eliminates the deep crevasses so the smoothed surface Hi,{x.) does 
have a skewed height distribution (with skewness dependent on the observation scale b). This is consistent with our 
own conclusion that the height fluctuations are non-Gaussian. 



B. Contour- loop analysis 



We perform loop measurements and check whether the different scaling relations derived in section IV are satisfied. 
The moments of the scale dependent curvature are used as an independent measurement of the roughness and to 
assess the Gaussianness of the height fluctuations. 

Using the loop algorithm (appendix ^), we measure the loop radii and corresponding loop lengths for 1000 contour 
loops constructed through randomly chosen points on the surface; however, we did not compute the loop correlation 
function. These loop measures support the scenario that the surface is self-affine up to a correlation length ^ w 25. 

The average loop length is plotted against the loop radius in Fig. |l^. We see a decade of power law scaling of the 
length with the radius, and from a linear least squares fit of the data to a line, for 5 < i? < 50, we find 

Df = 1.06(2) (7.2) 



for the fractal dimension of contour loops. Using the formula for Dj, Eq. ( 1.13| ), we calculate a — 0.88(4) in 



good agreement with the reported value, Eq. (^~^). In Fig. the dashed line corresponds to the percolation value 
Dh = 1.75; we see that loops at scales much larger then ^ show scaling consistent with this value. 

Finally, the number of loops whose length exceeds s, P>(s), is plotted in Fig. |l^. The data roughly shows two 
scaling regimes with different exponents, before it is cut-off by the system size. The knee occurs at loop lengths s ~ 70 
which, from Fig. |l^, corresponds to a loop radius of 20 or so; this again is comparable to the length scale ^ « 25 
found from the curvature data. From loops whose length is in the interval (10, 30) we extract the exponent 

T- 2 = 0.069(5) (7.3) 



while larger loops exhibit scaling consistent with the percolation value (indicated by the dashed line). Using Eq. ( 4.13| ) 
we find a = 0.85(1) again in good agreement with the self-affine exponent reported by Palasantzas and Krim. 

To summarize, the two loop measures we evaluated, as well as quadratic scale-dependent curvature, all indicate 
self-affine scaling with a roughness exponent a « 0.85, up to a length scale ^ « 25. Beyond this scale the height 
fluctuations appear to be uncorrelated. 



C. Other data sets 



The large value of a found for the silver-on-quartz STM data of Palasantzas and Krim is indicative of MBE-type 
growth which has surface diffusion as a dominant relaxation process. This motivates the study of other data sets 
which might correspond to different universality classes of growth. For example, the KPZ equation describes growth 
dominated by desorption and/or vacancy formation, both of which are relaxation processes that do not conserve 
particle number | |lO| , |ll|| . 

Loop measurements were carried out previously on gold electro-deposits by Gomez- Rodriguez et al pl[ |. These 
authors suggested the fractal dimension of contour loops as a useful measure for characterizing the surface morphology. 
What was lacking in their analysis was an equation relating Df to the roughness exponent a. From STM images 
of deposits grown in the fast and slow regime they determined the fractal dimension to be Df « 1.5 and Df k, 1.3, 
respectively. Now using Eq. ( 4.13 ) we calculate the roughness in these two regimes to be: a ~ and a w 0.4. The 
first is expected for Edwards- Wilkinson type of growth (a = 0), while the second is in good agreement with the 
Kardar-Parisi-Zhang value a = 0.38 (from most fits) or 0.40 (possibly exact) 47 4^. 

Csahok et al. studied the surface morphology of Ni films vapor-depositcd on a quartz substrate. (They were 
interested mainly in the effects of subsequent ion sputtering on the film.) We obtained an STM image of the as-grown 
Ni surface (before any sputtering) in the form of a 256 x 256 height array, and computed some of the contour-loop 
measures for it from a collection of 10000 loops. The results are consistent with a (KPZ-like) self-affine morphology 
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with a roughness of a = 0.4. Namely, the loop data shows a limited range of scaling for loop radii 10 < i? < 30. 
Direct fits to a straight line of the log-log plots of {s){R) and P>(s) in t he sc aling regime, yield Df = 1.31(4) and 
T — 2 = 0.22(2). After inverting the formulas for Df{a) and T{a) in Eq. (4.13) we obtain the estimates a = 0.38(8) 
and a = 0.44(4) respectively. 



VIII. DISCUSSION 



Here we summarize our main results, compare and critique previously introduced measures of surface roughness, 
and describe some open problems and interesting new directions. 



A. Summary of results 



We introduced (in Sec. Ill) new measures for characterizing the spatial correlations of rough interfaces. Their 
common property is that they are not linearly related to the structure function of the height. First, we introduced 
the scale dependent curvature. Its third moment is an indicator of the skewness of the height distribution, and thus 
is a good criterion for whether or not a surface's height fluctuations are Gaussian. Our chief focus, though, was 
on the ensemble of contour loops as a novel means of characterizing surfaces. For a rough self-affine surface, the 
loop ensemble is critical and we introduced three kinds of geometrical exponent associated with it: xi for the loop 
correlation function (probability that two points are on the same contour loop), the fractal dimension of a contour 
loop, Df, and r as socia ted with the length distribution of loops. In particular, we conjectured a super-universal value 
2x1 — 1 (see Sec. IV C ) which has been confirmed so far numerically (e.g. in Table ^, but not analytically. The 
loop exponents satisfy scaling relations (derived in Sec. iV), and granting the conjecture, their values are completely 
determined by the afline (roughness) exponent a. 

Next, we showed how numerical values of the geometrical exponents can be extracted in practice from height 
data obtained from simulations or experiments. We first did this in Sec. ^ for artificial Gaussian surfaces (known 
analytically to be self-afhne) and in Sec. VI for configurations from simulations of the single-step (growth) model 
(believed to be self-affine); this served as a check to confirm the validity of our scaling relations. Then in Sec. VII 
we processed an experimental data set - an STM image of a growth roughened silver film - in the same fashion. 
The results here also confirmed the scaling relations which in this case adds to the evidence of the self-afRne nature 
of the height fluctuations. The third moment of the scale-dependent curvature confirmed that the height fiuctuations 
are non-Gaussian, while the contour-loop fractal dimension and size distribution indicated self-affine scaling with 
a w 0.85. 

Experimental data often exhibit self-affine scaling up to a correlation length ^(i). We argued (in Sec. IV E| ) that 
the loop exponents, for loops whose linear size exceeds the correlation length, are determined by exactly known 
percolation exponents. The crossover between the self-afRne a nd p ercolative regime was visible (with a consistent ^ 
value) in every kind of measure on the expe rime ntal data in Sec. VII - the same was true for Gaussian random surfaces 
with an artificial length scale cutoff (Sec. VD). The numerical values of the percolative exponents were confirmed 
from simulations of Gaussian surfaces with a white-noise spatial power spectrum. 

Our results (see Sec. 0) show that it is quite difficult to get correct results from loop measurements when a is near 
to 1. The reason, we believe, is that the crossover to asymptotic behavior occurs at very large loops; the inferred a 
is thus smaller than the real one. It has been observed ||5^ that even the height-height correlation function tends 
to yield a too small value of a as compared to the Fourier power spectrum, even though the two measures have, in 
principle, the same information. 



B. Comparisons of roughness measures 

Roughness has often been analyzed based on a single number, the overall variance of the height over the entire 
system. However, spatial correlations in height fluctuations are central to the development of self-affine or other 
interesting morphologies. Therefore, every form of roughness measure we discuss takes the form of a spatial spectrum, 
i.e., one measures an entire function whose argument has dimensions of length (called s, R, b, or 1/q). The variation of 
the roughness measure with its argument is related to the varying amount of interface fluctuations on the corresponding 
length scales. 

Some previous measures of the self-affine exponents were reviewed in Ref. jist . They systematically compared the 
different measures using artificially constructed realizations of h{x) (only in 1 -|- 1 dimensions), and concluded that 
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the single best measure of a is the Fouri er p ower spectrum, our Eq. (^^). (Oddly enough, Ref. did not include 
the height correlation function, our Eq. ( p.2[ ), in their selection of measures to compare.) 

Another approach is to measure /i(x) along a single line in the x plane, corresponding to a line scan by the STM 
This section through the surface may then be analyzed as if it were a 1 + 1-dimensional profile. Ref. jsj] 
evaluated the variance over an interval of length Lqj which should scale as Lq"' ^^'^ applied this experimentally to 
the heteroepitaxy of CuCl on the (111) surface of CaF2. 

The most useful roughness measures have been discussed and critiqued in the sections related to them; they fall 



into three categories and are summarized here in Table III 



1. Quadratic roughness measures 



The most familiar measures are quadratic of which the first three were summarized in Sec. IIB 1 . Besides three 
well-known quadratic measures, we include a fourth which has not been previously applied: the variance of scale- 
dependent curvature, (Cfc(x)^), which we introduced in Eq. (|3.3| ) (Of course, the ensemble expectation cannot depend 
on x.) Its behavior is very similar to that of the height-difference correlation I?2(r), so (C^) is of interest mainly for 
comparison with the higher moments of Cf,(x). In practical applications the Fourier spectrum is probably the best of 
these. 

A key fact about the quadratic measures is that, given the complete function for any one of them, one can compute 
the complete function for any other one as a linear transform (convolution with some kernel) This property is not 
true for higher moments. Notice also that the quadratic measures are invariant under /i(x) —^ — /i(x) and so cannot 
possibly characterize the breaking of up/down symmetry in the growt h pro cess. Nor can they identify deviations from 



Gaussianness, since one can produce a Gaussian ensemble (as in Sec. VA) with any given Fourier spectrum 



2. Non-quadratic roughness measures 

Essentially all of these have been developed by analogy with quadratic measures, simply replacing the second 
power by a higher power. Our curvature-moment function seems to be the first generalization of the height-difference 
function that captures the up/down asymmetry. 

A simple generalization of the 6-box variance is the 6-box q-th moment, {{h{r) — hhY)b- The q = 'i moment 
characterizes the up/down asymmetry; when scaled by ((Wr) — hb)'^)^^^ it defines a scale-dependent, dimensionless 
skewness that measures the deviation from Gaussianness [pq ]. This appears to be a simple and attractive measure, 
but we know of no applications to date; our curvature moment (C^) is similar in spirit, but probably not linearly 
related. 

We evaluated the quartic curvature moment (C^), but this data was less useful than our other measures: it does not 
reveal the non-Gaussian nature as strikingly as (C^) does. The dimensionless ratio (C'^)/(C^)^ is 3 in the Gaussian 
case, but may not differ very much in a non-Gaussian ensemble. Furthermore, the roughness exponent was fitted less 
precisely from (C^) than from any other measure, probab ly due to the sensitivity of higher moments to rare events. 



The analysis in Ref. |pO| summarized in Sec. VII A 1 , appears to be the first application of a scale-dependent 



roughness measure to characterize the up/down asymmetry. However, we believe a local roughness measure such as 
the 6-box skewness or (better) our mean cubed curvature gives a more meaningful characterization. In a sense, the 
smoothed-height skewness is the opposite of the local measures since it includes the fluctuations from all length scales 
larger than b, while the local measures include the fluctuations from scales comparable to or smaller than 6; only the 
latter would be expected to scale as b". 

3. Loop measures 



The other non-quadratic roughness measures, of course, are the loop measures defined in Sec. III. The length and 
connectedness of a loop, manifestly, depend on the heights h{r) in a highly nonlinear fashion, and one might expect 
the loop exponents to be independent of the roughness exponent a; then the loop properties might have distinguished 
between different universality classes of growth which happen to have similar a values. From this viewpoint, it is 



disappointing that we in fact find the loop exponents are functions of a (Sec. IV). Thus for self-affine interfaces the 



loop measurements serve only as a check on other ways (quadratic and non-quadratic) of measuring a. Furthermore, 



when the heights at large separations are uncorrelated, implying the loops are percolation hulls (see Sec. IV E), the 
loop plots show a weaker change of slope at this crossover than the Fourier spectrum does. 
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It seems worthwhile nevertheless to compute loop measures. In a sense they depend on higher order correlation 
functions of the heights: then the agreement between the a values extracted from loops and from other measures is 
an additional, stringent test of self-afRneness. Also we observe empirically that loop measures, and in particular the 
average loop length as a function of loop radius, are very self averaging and measurement of a from them produces 
smaller errors than either the real-space or Fourier-space methods. Finally, although the loop exponents are the same 
for different universality classes with the same a, we do expect universal coefhcients to be different. 

For computer generated height data the loop-size-distribution is, perhaps, the single most valuable plot, because 
two different exponents can be obtained from scaling plots such as Fig. |^. This is not the case for experimental data 
where the system size is typically not a tunable parameter. 

The loop-correlation function G(y) is most tedious to compute, and since its exponent 2xi is superuniversal it does 
not yield an estimate of a. Nevertheless G(r) is a useful check on the self-afRneness, since the superuniversal behavior 
fails in other cases (e.g. beyond the correlation length, see Eq. ( f.l6| )). 



C. Future directions 



New experimental techniques which provide complete real-space images of the fluctuating quantity of interest (rather 
than system- wide averages, or local measures probing the system at only a few points), are being developed in every 
physical science. Consequently, measures which usefully exploit this wealth of information will gain in importance. In 
turn, the ability to measure new (and nonlinear) correlations may inspire new theories that can predict the correlation 
behavior. 



1 . Turbulence 



Fluid dynamics is a good example of the interplay just mentioned between theory and experiment: formerly two 
(sometimes more) point correlations were measured by hot-wire probes, and the same correlations were the objects of 
the Green's function method. As full images become available of the velocity field, many new measures are attempted 
in order to capture more of the available information. 

Indeed, the measures introduced here might be adapted to the geometrical description of turbulence. The advection 
of passive tracers by turbulent flows seems a to be an especially promising problem. There is already considerable 
interest in characterizing the equal-time correlations through fractal measures of the contours of (say) constant tracer 
concentration ||5^,^,^,||. (To maintain the analog of a surface's symmetry under global shifts of the height, one 
should study the logarithm of the concentration and use contours spaced equally on the logarithmic scale.) 

Measurements of the fractal dimension of iso-concentration lines of a passive tracer advected by a magnetically 
driven, turbulent, two-dimensional flow were reported by Cardosa et al. They find Df — 1.35(5), which, assuming 
the concentration field is self-affine, yields a roughness exponent a = 3 — 2Df — 0.3(1). Indeed, a = 0.30(3) was 
measured by the authors, by applying the q — 1 multiaffine correlation measure (entry 6. in Table |l|). We therefore 
infer that their measurements are consistent with a self-affine morphology for the concentration field. Details of the 
complete loop and curvature analysis of this data set will be reported in a separate paper |pq|. 



2. Other dimensions? 



In this connection, it is interesting to consider the generalizations of /i(x) to spatial dimensions of x other than 2. 
In the 1 -I- 1-dimensional case, there are no loops; the probabi lity o f first return to a fixed height value seems to 
be the closest analog to our loop size distribution of Subsec. [II C] (if presented as a distribution of R instead of s) 
and also to our loop correlation function Eq. ( 3.17 ). 

With each higher dimensionality there is greater richness in distinct geometrical measures that can be defined for 
iso-surfaces. For a hypersurface in 3+1 dimensions - like the concentration function in three-dimensional passive tracer 
advection - the level set may be multiply connected and even knotted. Nevertheless the size distribution exponent r, 
the fractal dimension Df, and xi of the connectedness correlations, can be generalized directly. But we see much less 
reason to expect a super- universal connectedness correlation exponent in dimensions higher than d = 2 + 1. 
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3. Multifractality and scaling relations 



Several mysteries reniain about the scaling relations der ived in Sec. IV. Above all, there is not yet any rigorous 
or analytic basis for our fundamental conjecture (Eq. ( 4.12| )) of a super-universal loop correlation that scales as 1/r 
for all rough self-affine surfaces - unaffected even by quenched disorder that further roughens the interface |Q. A 
second open question is to check numerically the correlation exponent a in Eq. (4.4) for an individual loop of fixed 
size; we did not evaluate it in any of our numerical studies, but it should not be the same as the exponent 2xi for 
the ensemble average over loops of all sizes. Finally, it is intriguing to ask what happens in a "multiaffine" system 
1^,^. Here different moments of th e hei ght variables have different scaling exponents; which of these (if any) is the 
one entering our formulas (4.9) and (4.10) for the loop exponents? 
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APPENDIX A: ANALYTIC DERIVATION OF XL{ct = 0) FOR AN EXACT SOLUBLE MODEL 



Our purpose here is to support the conjecture ( 4.12 ) in Sec. IV C , by showing that = 1/2 in the case of a lattice 
model which can be mapped to an equilibrium-rough surface. At long wavelengths height fluctuations are described 
by the well-known free energy 

F = (constant) y"d^x|V/i(x)|2 (Al) 

which by equipartition implies Eq. (|2.7|) with a = 0, so indeed the surface is self-affine. This appendix only summarizes 
arguments made previously in Refs. pl|| , , and . 

Consider a statistical model with microscopic heights zj defined on a triangular lattice {j}, such that zj changes by 
or ±1 between nearest neighboring sites. The partition function of the model is 

^=E n ^(^j-^k) (A2) 

{z} <j,k> 

where w(0) = 1 and w{±l) = K; the sum goes over all microscopic height configurations unrelated by a global height 
shift. 

A contour-loop configuration 7' is specified by drawing closed (periodic boundary conditions ensure that all contour 
lines are closed), oriented, non-intersecting loops along the bonds of the dual honeycomb lattice, which seperate sites 
that differ in height by ±1 (the sign determines the loops orientation). In terms of the loops, the partition function is 

Z = ^ X^" (A3) 

7' 

where K is the fugacity of an occupied bond (i.e. one covered by a loop), and Nb is the number of occupied bonds in 

. . . . . . . . 

This model is equivalent to the 0{2) loop model introduced by Nienhuis 133]. This is seen by rewriting the partition 

function in terms of non-oriented loop configurations 7, 

Z = Y,K'^'2^' (A4) 

7 

where Ni is the number of loops in 7 (which is the same as the number of loops in 7'), and the 2 appears as a result 
of summing over the two possible orientations for each loop in 7'. 
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By mapping the 0(2) loop model to the 4-state ferromagnetic Potts model on the triangular lattice, Nienhuis 
showed that for Kc = V2 the loop model is critical. Using the Coulomb-gas picture of correlations this implies 
that {zj} are rough; since the "background charge" is zero |Q, it is plausible that (at Kc) the height model is 



equilibrium-rough, i.e. the field /i(x) obtained by coarse-graining zj satisfies (Al). Furthermore the contour-loop 
correlation function can be identified in the 0(2) model with the so-called energy-energy correlation function which 
at K = Kc decays as a power law with the known exponent 

XI = 1/2 , (A5) 

as was to be shown. 

Since we view the 0{2) loop model simply as one of many possible lattice disretization of a logarithmically rough 
{a = 0) self-afhne surface, then the exponent xi = 1/2 should necessarily appear in other lattice models that map to 
rough surfaces. Indeed the same value of this exponent follows also from the exact solution of the 0{2) loop model 
on the square lattice and the n = 2 fully packed loop model on the honeycomb lattice 

APPENDIX B: PERCOLATION SCALING OF CONTOURS FOR UNCORRELATED HEIGHTS 

This appendix derives the scaling behavior of the loop ensemble when the random heights h(x.) have a finite variance 



and (beyond a correlation length ^(t)) are uncorrelated; this describes early stages of growth, as in Sec. IV E. 

To model the contour loops at length scales greater than £,{t), first coarse-grain the system into boxes of side ^(t). 
The average height h in each box is an independent random variable parametrized by p{h'), the probability that 
h < h'. 

Defining all the boxes with h < h' as "filled" simply reproduces the (uncorrelated) percolation clusters for occupancy 
p{h'). Then every contour of constant h' is simply the perimeter of such a cluster. This mapping is well-known from 
Rcf. and is widely applied in the theory of the quantum Hall effect [^,^ . 

The percolation clusters - as well as their perimeters - are self-similar only at p{h) — pc, the percolation threshold; 
we will first discuss their (known) loop exponents. The behavior when p ^ pc can easily be derived from well-known 
percolation scaling relations. The final step will be to integrate these results over p, since the loop ensemble we 
simulate actually corresponds to the union of perimeter ensembles for all p. 

1. Contour loops and critical percolation 

Fixing p = Pc for a moment, t he pe rimeter loop ensemble may be characterized by exponents Dh and Th, with 



definitions analogous to (3.9) and ( ^.13 ) for Df and r. (The subscript "/i" stands for "hull" as the perimeter is often 
called.) The fractal dimension 

Dh = 7/4 (Bl) 

is known exactly . 

The perimeter loops for percolation at pc also satisfy a hyperscaling relation analogous to (4.3), with a replaced 
by zero. That is, the largest cluster (or perimeter) diameter inside a box of side I is least ~ I. From this follows a 
relation for r/j 

Th = l + 2/Dh - 15/7 . (B2) 

When p ^ Pci the cluster (and perimeter) ensemble scaling is cut off at the percolation correlation length £,p{p), 
which diverges near pc as 

Up)-\p-pcr'' ^ (B3) 

where Vp is the usual percolation correlation exponent, and = 4/3 is known exactly |^,^^. In this case, the loop 
length distribution is 

P,is;p) = s-^^--^^Ms/ip{p)^-) (B4) 
where fh{) is a scaling function, which falls off exponentially fast for loops of radius greater than ^p{p). 
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2. Union of all percolation contours 



In the percolation regime, evidently, the statistical properties of the contours of a particular level set depend on the 
chosen level h. (This was impossible in the self-afRne regime, since in that case the fluctuations of h were unbounded.) 
But we have previously studied the union of all contours with different h, c orres ponding to all values of p{h) from 
to 1. That is, indeed, the ensemble sampled by our computer codes (s ee Se c. VC ). W e will now derive the exponents 
Tp and xi^p of this ensemble, defined analogously to r and xi in eqs. ( ^3.13 ) and ( 3.17j ). 

Most of the loops at a large length scale R come from levels sets at height h with £,{p{h)) > R, rather than from 
the exponential tails of the distribution (B4) for the other h values. Thus these obey the percolation scaling and all 
have fractal dimension Df^p = Dh- This behavior is illustrated in Fig. ||(c). 

Now, in the percolation regime, P{s) as defined in Sec. Ill C , is just proportional to the integral of Ph{s;p{h)) over 
h. A weighting factor \dp/dh\ should be included as the contours are equally spaced, and p{h) is normalized to unity. 
Since the large contours come from p ^ pc, only that part of the distribution matters. Inserting (BS) into (p34|) , one 
obtains 



P{s) ^ / dp s-'^^'^-'^ h{iconst)s\p ^ p.f'^-'-) 



(B5) 



hence P{s) ^ s '"^p with 



Tp = Th + {DhVp)-^ = 18/7 



(B6) 



Finall y, gi ven (B£) the simplest route to the loop (connectedness) correlation exponent is to use the exponent 
relation (|4.8|); this gives 



2Dh + j^: 



5/4 



(B7) 



Eq. (B7) could alternately be reached by first noting that the corresponding exponent is 1/2 for the percolation hull 
ensemble at Pc, and then averaging the loop connectedness correlation function analogous to (|B5|). 



APPENDIX C: LOOP FINDING ALGORITHM 



Given a square lattice C on which the heights h are defined and a point Xq on the dual lattice £*, the task of the 
loop finding algorithm is to construct a contour loop of the surface which passes through the point xq . 

The contour is a walk along the bonds of C* that cuts those bonds of C that have vertices with heights lying above 
and below the contour height. Fig. |l8|. To implement this idea we first define the level height /iicv which is the average 
of the four heights around the plaquette centered at Xq. Second, we assign to all the sites of £ + or — signs according 
to whether they are above or below the chosen level /iiov Now, starting from Xq we form the contour loop by drawing 
links on the dual lattice which cross the bonds of C connecting + and — sites. This is repeated until the walk returns 
to the starting point xq; the finite extent of the lattice C is dealt with by implementing periodic boundary conditions. 

Special care must be taken whenever a "saddle-point" plaquette is reached, that is one where the sites of the lattice 

are assigned H \ — signs cyclically around the plaquette. In this case four links meet at the point in the center and 

we must resolve the connectivity there by an additional rule, so as to convert this pattern into two 90-degree turns 
that are not quite touching. One natural condition on the rule is that it should be reversible, that is one should find 
the same loop whether one starts traversing it clockwise or counterclockwise. A second condition is that it ought to 
be invariant under reflecting all heights by /i(x) — > — /i(x). A physically sensible rule which satisfies both conditions 
makes use of the average height /ipiaq of the four heights around the saddle-point plaquette. If /ipiaq < /iiev, we view 
the center of the plaquette as being lower than the level of the contour loop and the connectivity is resolved by having 
the -|- sites inside the 90-degree turns. In the opposite case, /ipiaq > /iicvj the -I- sites lie outside the 90-degree turns; 
see Fig , [l^ . (The agreement of loop data from the single-step model with parameters p+ and with 1 — p+, as explained 
in Sec. VI was a valuable check of the up-down symmetry of our loop-finding algorithm.) 

Once a contour loop through Xq has been constructed its length and radius are recorded, assuming that the loop 
is topologically trivial. (Due to periodic boundary conditions loops with non-zero winding numbers are possible and 
these we discard.) The contour loop length s is equal to the number of steps made during the loop construction, 
while the radius R is the size of the largest square which covers the loop, i.e., the maximum displacement in the x 
or y direction. Every topologically trivial loop also contributes to the correlation function G'(r); for every point on 
the loop that is a distance r € [i,i + \) {i is an integer) away from the starting point xq, the array element g{i) is 
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increased by one. We def ine G{r) for our simulations as g{r)/2T:r, which asymptotically is normalized the same as 
G(r) defined in Sec. |lliq . 

So far we have assumed that the height variables are real and the condition /ipiaq = /iicv is almost never fulfilled. 
This is not the case for interfaces which arise from discrete growth simulations like the one presented in section ^ 
where the height variable takes on integer values. In this case the resolution of the connectivity should be completely 
random but we must ensure that we use the same choice if the loop returns to the same plaquette. The simplest 
way to do this, which is what we have implemented, is to take the original integer heights and "dither" them - add 
small amounts of random, uncorrelated gaussian noise to all /i(x). This will also solve the problem of choosing /iicv! 
it will be non-generic for any two heights to precisely coincide, although this will happen occasionally as the price of 
roundoff error. When this does happen we start over by choosing a new initial site Xq. 
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TABLE I. Geometrical exponents xi, Df, and t for loops on Gaussian surfaces with various roughness exponents a. Columns 
marked "direct" are from direct fits to a power law of the data from system size L — 512, inferring '2x1, Df, and r — 2 from plots 
such as Figs. ^, and Columns marked "FSS" were fitted to finite-size scaling plots like Figs. g(b) and |^(b). According 
to our conjecture, the "theory" value of 2xi is 1, independent of a, and this is supported by the measurements here. Notice a 
slight systematic deviation of the "direct" exponents from theory when a > 0. 5, which we attribute to more severe finite size 
effects in those cases. The "theory" formulas for Df and r — 2 are in Eq. (4.13). 
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2x1 


Df 


r - 2 1 




direct 


FSS 


direct 


FSS 


theory 




direct 


FSS 


theory 


0.0 


1.07(2) 


1.02(2) 


1.48(1) 


1.50(2) 


1.5 




0.35(2) 


0.33(1) 


0.333 . . . 


0.2 


1.04(1) 


0.97(2) 


1.39(1) 


1.41(2) 


1.4 




0.30(1) 


0.28(1) 


0.286... 


0.4 


1.01(1) 


0.98(2) 


1.31(2) 


1.32(3) 


1.3 




0.24(1) 


0.225(5) 


0.231... 


0.6 


1.00(1) 


0.97(2) 


1.23(3) 


1.19(3) 


1.2 




0.18(1) 


0.165(5) 


0.166... 


0.8 


0.97(2) 


0.97(2) 


1.15(1) 


1.11(2) 


1.1 




0.12(1) 


0.11(2) 


0.090 . . . 


1.0 


0.95(1) 


0.96(2) 


1.06(2) 


1.04(3) 


1.0 




0.08(2) 


0.02(2) 


0.000 . . . 



TABLE IL Results of fits to roughness measures applied to surfaces generated by the single-step model. The exponent values 
a in the left two columns were derived in two direct ways, from the data plotted in Figures |l|(a) and The direct-fit results 
("direct") for Df and r — 2 used only the L = 128 data (shown in Figures |l^ and ^J); the finite-size scaling results ("FSS") 
were obtained using system sizes L = 32,64, 128. The su bhea dings "a" under Df and r — 2 are estimates of the roughness 
exponent obtained from the FSS results by inverting Eq. (4.13). 
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direct 


FSS 
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direct 


FSS 
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0.1 


0.33(2) 


0.35(1) 


1.38(1) 


1.35(2) 


0.30(4) 


0.30(1) 


0.24(1) 


0.37(3) 


0.3 


0.19(4) 


0.09 (2) 


1.47(1) 


1.46(2) 


0.08(4) 


0.38(1) 


0.35(1) 


-0.08(5) 


0.5 


0.135(2) 


0.08 (1) 


1.51(2) 


1.50(2) 


0.00(4) 


0.40(2) 


0.36(2) 


-0.13(5) 



TABLE in. Roughness measures 



Quantity 


Description 


Quadratic measures 


1. {(ft(x)-M')6 


variance in 6 x 6 patch 


2. Z?2(r) = (|Mr)-M0)|^) 


height correlation 


3. {|Mq)P) 


Fourier spectrum 


4. (la(x)i^) 


"&-dependent curvature" variance 


Cubic and other non-quadratic measures 


5. {{h{r) - htf)t 


skew moment in & x b patch 


6. {\{h{r) - h{0)n 


g-multiaffine height correlation 


7. {{H.ir) - hf) 


skew moment of 6-smoothed height 


8. {[a(x)]^) 


"curvature" skew moment 


9. {[^(x)]*) 


"curvature" quartic moment 


Loop measures 


10. {s)r 


average loop length, given radius R 


11. P>(s) 


Prob (loop through x is longer than s) 


12. G(r) 


loop connectedness correlation function 
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FIG. 1. Contour plot of a a = 0.4 random Gaussian sur- 
face; system size L=512. 



26 




(s> 



PJs) 



-.1 ■75/ 



0.57 




FIG. 2. Percolation crossover. In each figure, tlie solid line 
represents the function at a certain time t, and the dashed 
line represented the same function at a later time when ^{t) 
has increased. The exponents shown in the figures are for 
a = 1/2, but the qualitative behavior is the same so long 
as < Q < 1. Each graph shows a crossover to percolation 
exponents at a "knee" which corresponds to a length scale 
~ £,{t): (a) Fourier spectrum (|ft(q)|^), (b) loop correlation 
function G(r), (c) average loop length < s > (R), and (d) 
cumulative distribution P>{s) of loop lengths (through a given 
point). 
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FIG. 3. Scale-dependent curvature moments from Gaus- 
sian random surfaces with roughness exponents a = 0.4 (cir- 
cles) and 0.8 (triangles). The third moments (lower plot) 
are zero confirming the up/down symmetry. The upper plot 
shows the second moments (open symbols) and fourth mo- 
ments (filled symbols). 
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FIG. 4. Average loop length (s) as a function of loop 
radius R, for random Gaussian surfaces with a = 0, 0.4, 0.8 
(from top to bottom); system size L = 512, and 10* loops 
were collected. The "direct" Df data in Tableware obtained 
by linear least-squares fits to such plots in the scaling regime, 
which is roughly 10 < i? < 100. 




FIG. 5. Cumulative number of loops whose length is big- 
ger than s for random Gaussian surfaces with a = 0, 0.4, 0.8 
(from top to bottom); system size L — 512. Here and in all 
other plots of P>(s), raw data is binned in intervals of form 
(s,l.ls). 



28 




FIG. 6. Loop correlation function for random Gaussian 
surfaces with a = 0,0.4,0.8 (from bottom to top); system 
size L = 512. In this and all such plots, raw data is binned 
logarithmically in intervals of form (r, l.lr). The latter two 
graphs are offset vertically by factors of 10 for clarity; they 
arc virtually identical except for a "knee" at slightly different 
r values. 
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FIG. 7. Cumulative loop-size distribution, for random 

Gaussian surfaces with q — 0.4. (a) Data for system sizes 
L = 64 (O), L = 128 (A) L = 256 (v), and L = 512 
(*). (b) Collapse of this data in a finite-size scaling plot with 
r - 2 = 0.225 and = 1.32. 
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FIG. 8. System-size dependence of the loop correlation 
function G{r) for Gaussian random surfaces with a = 0.4. (a) 
Data for sizes L = 64 (Q), ^ = 128 (v) L = 256 (v), and 
L = 512 (*). (b) Data collapse of this data with 2xi = 1.02. 
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FIG. 9. Fractal dimension Df and length distribution ex- 
ponent r — 2, as functions of the roughness exponent q of a 
random Gaussian surface, obtained from finite-size scaling fits 
(see Table ^. T he so lid and dashed lines corresponds to the 
formulas in Eq. (4.13). 
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FIG. 10. Finite correlation length effects for Gaussian ran- 
dom surfaces with a — 0.4, and a crossover to white noise for 
wavevectors smaller than ir/^q, where = 16; L — 512 is 
the system size. Circles are used for data sets with no cutoff 
which are included here for comparison with the cutoff data 
(triangles), a) Squared curvature function - note the knee at 
b ~ 15. b) Average loop length as function of radius - knee 
at ii ~ 20. c) Cumulative distribution of loop sizes - knee at 
s ~ 100. d) Loop correlation function - knee at r « 20. 
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FIG. 11. Power spectrum of the height in the SSM model 
along the [1,0] (filled symbols) and the [1, 1] direction (open 
symbols) in reciprocal space. The data for p = 0.3 (triangles) 
and for p+ =0.5 (diamonds) has been shifted with respect to 
the p+ = 0.1 data (circles) by factors of 0.1 and 0.01 respec- 
tively (for clarity) . Note that the power spectrum is isotropic 
in Fourier space for small values of |q|. 
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FIG. 12. Second moment (a) and third moment (b) of 

the scalc-dopcndent curvature C6(x), for surfaces from the 
single-step model with p+ = 0.1 (A), p+ = 0.3 (o), and 
p+ = 0.5 (*). In (b), the p+ = 0.3 and 0.5 data are con- 
sistent with (Cj) = 0, while the = 0.1 data show a strong 
(and non-Gaussian) breaking of up/down symmetry. 
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FIG. 13. Average loop length "vs" the loop radius for 
p+ = 0.1 (A), p+ = 0.3 (o), and p+ = 0.5 (*), in the sin- 
gle-step model. Note that the p+ = 0.3 and p+ = 0.5 data 
are almost indistinguishable. 
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FIG. 14. Normahzcd loop size distribution for p+ = 0.1 
(A), p+ = 0.3 (o), and p+ = 0.5 (*), in the single step model. 
Again, the p+ = 0.3 and p+ = 0.5 plots are almost indistin- 
guishable. 
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FIG. 15. Second moment (a) and third moment (b) of the 
scale-dependent curvature, as evaluated for a 702 nm thick 
Ag film, grown on quartz, from the STM data of Palasantzas 
and Krim (Ref. M). 
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FIG. 16. Mean contour length (s) as a lunction of radius 
R, for the Ag film of Ref. ^9|. Here 1000 contour loops were 
coUected from the STM data of Ref. @. The solid hne is the 
least-squares best fit for radii 2 < R < 125; its slope is the 
estimated fractal dimension Df. The slope of the dashed line 
is equal to the hull dimension of critical percolation clusters. 
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FIG. 17. Cumulative distribution ol contour loop lengths 
from STM data, for the Ag film of Ref. jig). The solid line is 
the result of a linear fit to the data in the affine-scaling regime. 
The slope of the dashed line corresponds to the exponent r — 2 
in the percolation regime. 
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FIG. 18. Construction of contour loops of a random sur- 
face on a lattice. Heights h{r) are indicated by numbers in 
the cells; /iiev is the height of the level set through the cho- 
sen point (filled circle) while /ipiaq is the height of the "sad- 
dle-point" (unfilled circle). Our definition of the diameter 
R and loop length s is indicated. The solid arrow connects 
points on the same loop, and thus contributes to the loop cor- 
relation function G(r); the dashed arrow does not contribute: 
it connects points of the same level set, but they are on dis- 
connected loops. 
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